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Abstract 

Exact and closed-form expressions of the particle density, the kinetic energy density, the proba- 
bility current density, and the momentum distribution are derived for a coherent state of a nonin- 
teracting Fermi gas, while such a state can be obtained from the ground state in a d-dimensional 
isotropic harmonic trap by modulating the trap frequency and shifting the trap center. Conserva- 
tion laws for the relations of the densities are also given. The profile of the momentum distribution 
turns out to be identical in shape with that of the particle density, however, the dispersion of 
the distribution increases (decreases) when that of the particle density is decreased (increased). 
The expressions are also applicable for a sudden and total opening of the trap, and it is shown 
that, after the opening, the gas has a stationary momentum distribution whose dispersion could 
be arbitrarily large or small. 
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Experiments on quantum degenerate Fermi atoms provide motivation for the theoretical 
analysis of the noninteracting Fermi gas in a harmonic trap. Since the interaction energy 
per particle is small compared with the Fermi energy, the interactions between atoms are 
controllable, and since the short-range interaction may be ignored in the long-term behavior 
of the expanding Fermi gas, theoretical studies of noninteracting particles are important for 
understanding experiments on degenerate Fermi gases Q|. Indeed, in the deep BCS limit, 
the Fermi gases of atoms are considered to be noninteracting 0, y]- As most experiments 
are performed in harmonic traps, analytic expressions of the particle density and the kinetic 
energy density have been given for the (M + 1)— shell filled ground state of the Fermi gas 
in a d-dimensional isotropic harmonic trap P|. 

For an isotropic harmonic oscillator, the time-dependent Schrodinger equation is invari- 
ant under the displacement-type transformation (DTT), and under squeeze-type transfor- 
mation (STT) with a rescaling of time. As implied by Kohn's theorem, the system of a 
single-component gas is still invariant under the extended DTT, which can be used to find 
a coherent state of the harmonic motion of the center of mass p, Iq. The system of a non- 
interacting gas is invariant under the extended STT, which predicts a generalized coherent 
(squeezed) state of breathing (compressional) motion 0, Q, • In experiments, the system 
of the (generalized) coherent state can be prepared by exciting the system in an eigenstate 
through shifting the trap center and modulating the trap frequency. For the one-dimensional 
boson gas of infinitely strong repulsive interaction of zero-range (Tonks-Girardeau gas), the 
particle density is identical to that of a noninteracting Fermi gas, and the invariance un- 
der the STT may be applicable to some extent 0, 1^ . Recently invariance under the STT 
has been used to analytically explain the "dynamical fermionization" of the momentum 
distribution of the expanding Tonks-Girardeau gas [gl llo|. 

In this paper, exact and closed-form expressions for the particle density, the kinetic energy 
density, the probability current density, and the momentum distribution are presented for a 
state coherently excited from the (M + 1)— shell filled ground state of a noninteracting Fermi 
gas in a (i-dimensional isotropic harmonic trap, and conservation laws for the relations 
of the densities are given. The particle density and the momentum distribution, which 
have been studied as a function of the scattering length through a mean-field theory in 
Ref. j3[, are of great experimental interest , and the expressions presented here are 
also valid for the analysis of a sudden and total opening of the trap. The profile of the 
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momentum distribution turns out to be identical in shape with that of the particle density, 
however, as a manifestation of the Heisenberg uncertainty principle, the dispersion of the 
distribution increases (decreases) when that of the particle density is decreased (increased). 
The dispersion of the distribution could be arbitrarily large or small, as the amplitude of 
the breathing motion could be jQ]. Though only the coherent states excited from a ground 
state will be considered here, since the large or small dispersion basically comes from the 
invariance under the STT and the uncertainty principle, general coherent states could also 
have the large or small dispersion. After a total opening of the trap, we show that the 
momentum distribution of the free gas remains stationary without changing from that of 
the gas at the moment of the opening. Since we ignore the interactions, we consider only 
the spin-polarized case. 

For a d- dimensional isotropic harmonic oscillator described by the Hamiltonian 

Hiv, p, t) = ^ ( |- + -mw'{t)x^ - F,{t)xi ) (1) 

with m, w{t), and -Fj(t) denoting the positive mass, real frequency, and an external force 
in the i-th. direction, respectively, through the extension of the one-dimension result 
a complete set of the wave functions ipn{j,t) satisfying the time- dependent Schrodinger 
equation is found, with n = (ni, n2, . . . , n^), as 

d 

Mr,t) = l[^nXXi,xlt), (2) 
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where rii is a non-negative integer, and 
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with u, V being two linearly independent real homogeneous solutions of the classical equation 
of motion 



Xcl + W^(t)Xcl 



m 



and ri{t), I being defined as 



(4) 
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respectively, while a positive constant Q is defined as = m[vu — iiv], and the overdots 
denote differentiations with respect to t. In Eq. (3), is the Hermite polynomial, 
denotes a solution of Eq. (4), and 6i is defined through 6i = y [w^(t)x^ — x^] . 

When w(t) is a constant Wc and Fi{t) = 0, if we choose u = cos{wct), v = sm{wct), 
and xf = 0, then / = a/ h/ (mwc), and tpni^^t) is the wave function of an eigenstate. The 
(M + 1)— shell filled coherent state can be obtained from the (M + 1)— shell filled ground 
state, by applying an external force Fi{t) and modulating w{t). If the coherent state is 
established once, the system remains a coherent state in general even after the external 

one-particle reduced density 

matrix of the (M + 1)— shell filled coherent state is given by 



p(r',r;t)= C(r', t)V'n(r, t). 

With an independent variable b, we define 

{n} 

Making use of this fact 
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which comes from Mehler's formula 
C(r',r;t;6) 



^2 sinh(6) 



12 1, we find that 



(6) 



(7) 



(8) 



where 



with 



X exp 
X exp 




Fp ■ S + 



tanh I - 



77?coth 
4/2 



q 



r)/2, 



p 



15 -^25 ■ ■ ■ ) 



Just as in Ref. P], the infinite summation in Eq. ((Tj) can be 



(9) 
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truncated, through the well-known facts in the (inverse) Laplace transformation 



A [e{p-k)] = e-'''/b 



with a positive constant k and the step function 0(s) satisfying 6(s) = 1 for s > and 
G(s) = for s < 0. The density matrix may thus be formally written as 



p(r',r;t) = /:-i 



^C(r',r;t;6) 



(11) 



with X = M + {d + l)/2. By explicitly carrying out the inverse Laplace transformation for 
the diagonal component of the density matrix, we find the particle density 

pir,t) = pir,r;t) = ^Viz), (12) 

where 

z = {r-r,r/f (13) 
and the function 'D{z) is defined as jj] 

^(^) = ^E(-l)"^M-n^"W' (14) 
n=0 

with = {d/2 + 2n)T{d/2 + n)/[n\{d/2)\], F^^^ = 2T{d/2 + n + l)/[n\{d/2)\] and L^{z) 
being the Laguerre polynomial. 

The probability current density defined by 

S(r, t) = (^n(r, t) VV^n(r, t) - ^„(r, t) VC(r, t)) (15) 

with denoting the truncated summation of Eq. ® , can be found as 

S(r,t) = Av,p(r',r;t)|,^o = vp(r,t), (16) 
im 

where 

v = r,+ (/7/)(r-r,). (17) 

Up to the phase, the wave function ■?/'n(r, t) can be obtained from that of an eigenstate by 
isotropically rescaling the space according to the ratio I and globally displacing it by the 
amount r^. Since, under the rescaling, /// plays the same role as the Hubble constant in 
the expanding universe, v is the velocity of the position r under the transformations, which 
explains the classical aspect of S(r,t). 

For the kinetic energy density, we explore the quantity 

T-s(r,t) 
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-2(VC(r,t))-V^„(r,t)], 



(18) 



which can be evaluated as 
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P(r,t) + rQ(z,t), 
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where 



coth(|) exp[— 2;tanh(|)] 



By carrying out the inverse Laplace transformation explicitly, we find that 



rQ{z,t) 



frde 



[Z 



n=0 



with G^l = (32ri2 + (P + 16nd + 2(^)r(d/2 + n)/[4{n\)T{d/2 + 2)], and G^^^+i = 2(4n + ci + 



(20) 



(21) 
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2)T{d/2 + n + l)/[n\T{d/2 + 2)] TQ{z,t) can be represented in terms of p(r, t) as 



'rQ{z,t) = 
d + 2m 



_VV(r,t) + ^ ^ -Pir,t) 



(22) 



which, for the ground state, reduces to the known relation An expression of the 

kinetic energy density is then written as 



2m; 

mv^ , , ^d + 1 , , 
^P(r,t) + 2^rQ(z,t) 



p(r', r; t) satisfies the following Schrodinger equation 
d 



= -|^(V^-V'^)p(r',r;t) 



+ 
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r — r 



[r - r') ■ F(t) 



p(r',r;t). 



(23) 



(24) 



6 



which can be used to find conservation laws [1J| . The coincidence hmit r — r' of Eq. 
gives 

^p(r,t) + V-b(r,t)v] = 0. (25) 

The coincidence hmit after applying Vs on Eq. yields 

d d Ax' d 

= -p{r^t)[w\t)x,-^]. (26) 
Similarly, by applying Vg on Eq. with the help of Eq. (fTT|) , we find that 

df"^^^'^^ + [^q(^>^)v.] = -2jrQ(2;,t). (27) 
^^'or the ground state of w{t) = Wc and Fi(t) = 0, Eq. reduces to the known relation 

The momentum distribution defined by 

n{p, t) = ^ I drdr'e-'^<^-^'^/''p{r', r; t) (28) 



can be calculated as 



(27r)'^/2 



-[(p-mrp)2tanh{|)/a2(t)] 



6sinh'^/2(6) 

= a~\t)V{{p-mi,f/a\t)), (29) 

where 

a(t) = ^hyP + m'^P. (30) 

The momentum distribution in the momentum space thus has the identical shape of the 
particle density in coordinate space upon di^lacement and rescaling, as has been estimated 
for the one-dimensional free expanding gas [8]. Just as in the particle density (see, e.g., 
IigI). with the particle number N, the momentum distribution goes over to the approximate 
form 

_ [2{d\NY/'-a-\t)ip-mr,yf' 
"^TF 2^7r'^/2r(l + d/2)a'^{t) ' ^ ' 

for 2{d\NY^'^' > (p — mrp)2/a2(t), in the larg^iV limit, and the exact profile has corrugation- 
like corrections from the approximation [lly]. For a time-constant trap, if we choose u = 
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cos{wct), V = sm{wct), and xf = 0, utf exactly agrees with the distribution found for the 
deep BCS hmit of the atomic Fermi gas through the mean-field theory ^ (see also Ref. 

For a general w{t), I satisfies ml + mw'^{t)l — f}} /{ml^) = 0, and the behavior of / has 
been analyzed for time-periodic wit) [ll| to find that I could be more and more amplified 
in general and become unstable as time passes. For a time-constant trap, as in the case 
after the modulation of the frequency is stopped, / satisfies m^/^ + h^/P = —rn?w1l'^ + 
2me = a^(t) with a constant e, which shows, as a manifestation of the uncertainty principle 



17[, a(t), characterizing the dispersion of the momentum distribution, increases (decreases) 
when / or the dispersion of the particle density is decreased (increased). For the time- 
constant trap, up to a global time-displacement, / is written without losing generality as 



/ = a/V {mwc)^/ A? cos^ w,t + sin^ WftjA^^ with a positive constant A [6j. In this case, 
a(t) = a/ mhWf\A^ sin^ w^t + cos^ Wc^/^^], which shows, the dispersion of the momentum 
distribution could be arbitrarily large or small in the limit of A — or A — oo. 

The formalisms developed are also valid for w{t) = in d-dimension, with u = ai, 
V = a2{t — to) + as of 01,02,03,^0 being constants. For a sudden and total opening of the 
trap at t = to, the continuity of the wave functions requires the continuities of , xf , rj, rj, 



which gives 02 = a/^^^/ {mr]_.Y + ^-5 '^1 = ^7^5 '^3 = , where the subscript — denotes 
that the quantity is evaluated in the limit of t — to ~^ ~0. If ?)_ < 0, the dispersion of 
the particle density, which is proportional to Z, will decrease after the opening of the trap 
until t — to = |a3|/a2, and then it will increase, while the cloud of the gas will continuously 



expand for ?)- > 0. a{t) of Eq. (p?n|) with / = ^/h[af + (a2(t — to) + cl^)"^] / {maia2) is a 



constant a;_, that is, for t > to, a{t) = ^Jmha2/ai = yh^/l"^ -\- m'^P_ = a-. After the 
opening, the momentum distribution of the gas is thus stationary, as has been predicted for 
the one-dimensional expanding gas. If the opening is for the ground state of the trap of a 
constant frequency Wc, I- = \/hj (mwc), and /_ = 0, and thus a{t) = y/hmwc ^. 

In summary, the expressions of the particle, kinetic energy, probability current densities, 
and the momentum distribution have been derived for a coherent state of a noninteracting 
Fermi gas in a (i-dimensional isotropic harmonic trap, and, for the relations of the densities, 
conservation laws have been given. It has been shown that the profile of the momentum 
distribution is identical in shape with that of the particle density, while the dispersion of the 
distribution increases (decreases) when that of the particle density is decreased (increased) in 
a time-constant trap. The dispersion of the distribution could be arbitrarily large or small, 
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which basically comes from the invariance under the STT and the uncertainty principle. 
Since the invariance exists only for restricted cases , an appearance of the distribution 
of the large or small dispersion in a system of Fermi atoms may indicate that the system is 
in the noninteracting regime. After a total opening of the trap, the momentum distribution 
does not change from that of the gas at the moment of the opening, while, long after 
the opening, the dispersion of the particle density increases according to a^t/m with the 
constant a;_ which characterizes the dispersion of the distribution. As the uncertainties are 
controlled by the classical solutions in the coherent states, an experimental realization of 
the system described here may be important to provide an observable manifestation of the 
uncertainty principle. 
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